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Abstract 

In some respects, a cluster consisting of many atoms may be regarded as a single large atom. 
Knowing the dielectric properties of such a cluster permits one to evaluate the form of the van der 
Waals (dispersion) interactions between two clusters or between one cluster and a surface. In this 
paper, we derive these interactions in two extreme opposite regimes of separation: fully retarded and 
nonretarded. In the fully retarded regime (very large separation), the magnitude of the interaction 
is determined by just the static polarizability of the cluster(s). In the nonretarded regime (small 
separation), we employ a single resonant frequency model of the cluster polarizability to derive 
expressions for the interactions' coefficients. Numerical examples are presented to demonstrate that 
many-body screening of these interactions can be significant. The results represent the corrections 
to the commonly used approximation of pairwise additivity of interatomic interactions. 
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I. INTRODUCTION 



Van der Waals (vdw) interactions have been much studied in recent years, in part be- 
cause of increasing interest in phenomena involving surfaces, interfaces, and nanoparticle 
interact ionsSSffli'lfl. In this paper, we evaluate two distinct interactions involving clusters, 
which we assume to be spherical. One of these is the interaction between two clusters, de- 
noted A and B; the other is the interaction between a cluster and a flat surface. We consider 
both of these interactions in two extreme opposite regimes of separation: nonretarded and 
fully retarded, corresponding to small separation and very large separation, respectively. The 
word "small" is used advisedly because we assume throughout this paper that the separation 
is large compared to the size of the clusters. 

Nonretarded vdw interactions, sometimes called London forces, can be derived from quan- 
tum mechanics and electrodynamics. Most research involving interactions within condensed 
states of matter pertains to the nonretarded regime because short-range forces predominate 
in determining the physical behavior of such systems. However, the onset of effects of re- 
tardation occurs at sufficiently small distance scales (~ 10 nm) that such effects are often 
present to some extent; they become particularly important when studying large molecules 
or clusters, because of their size. The interaction is said to be fully retanie<iH!0 in the 
separation regime where the finite speed of light causes a significant attenuation of the 
(nonretarded) vdw interaction. Indeed, even the power law describing the distance depen- 
dence of the interaction is altered by retardation (increasing the falloff by one power of 
distance). The regime of retardation is that where the separation r is large compared to 
both the cluster size (radius R) and a typical wavelength A determined by the spectra of 
both interacting species: A = he/ AE; here AE is a characteristic energy in the electronic 
spectrum of the cluster and/or surface. If AE is the first ionization energy, for example, 
then A ~ 0.1 to 1 micron (for Na and Si, A ~ 0.24 microns and 0.15 microns, respectively). 
The surface force apparatus has been used to measure van der Waals forces between two 
macroscopic mica surfaces^. Other techniques (e.g., atomic force microscopy^! and total 
internal reflection microscopy^! have been used to measure interaction forces, and a few tech- 
niques (e.g., laser trapping^! and differential electrophoresis!)!) have been used to measure 
sub-piconewton interaction forces between Brownian particles. Figure 1 presents results for 
the interaction between two Na clusters obtained in this paper for the nonretarded and fully 
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retarded regimes, along with a plausible interpolation over the interval 0.1 < r/A < 10. 

Because vdw interactions arise from fluctuating electromagnetic fields in the interacting 
media, there arise many-body effects associated with screening of these interactions by 
neighboring particles within a cluster. This means that one should not simply add together 
the interactions of individual particles without taking the screening into account@. As an 
example, we demonstrate the importance of screening in the case of alkali metal clusters. 
There occurs a general increase in the effect of screening as a function of cluster size, but 
the trend is not completely monotonic because of electronic shell effects on the cluster's 
polarizability. 

The outline of this paper is the following. Section II describes the case of fully retarded 
interactions. In that case, the vdw interaction is determined by just the static polarizability 
of the interacting clusters (because long wavelength fields give rise to the interaction). In 
Section III, we explore the nonretarded regime, for which all wavelengths contribute, in 
principle. That treatment requires the introduction of a model form for the polarizability 
as a function of frequency. In the present case, we employ the simplest possible model 
(single resonant frequency) that is consistent with both the static value and the known high 
frequency limit of the polarizability. Section IV summarizes our results and comments about 
possible extensions to related problems. 



II. FULLY RETARDED REGIME 



The starting point of our analysis is the fact that the fully retarded interaction between 
two atoms (here denoted a and b) separated by distance r (with r ^> A) is known to satisfy!: 

P°P° 

V ab (r) = -K^ , K = ™£ (1) 
r' oir z 

Note the interesting fact that the interaction falls off with distance as r -7 , with a co- 
efficient determined by the static atomic polarizabilities (P°, P b °) of the species involved. 
This behavior contrasts with the nonretarded vdw interaction in both the latter's power law 
(r -6 ) and its dependence on the frequency- dependent polarizability, rather than just the 
static value P° (see Eqs. |1(] and [11] below). 
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A. Cluster-cluster interaction 



We now turn to the problem of the interaction between two clusters. The key idea of this 
paper is simple: one may treat a cluster as though it is a large atom, insofar as the cluster 
has an excitation spectrum with an energy gap and it is of finite extent. Here, we consider 
the case of separation r large compared to its size. Specifically, the present calculation 
requires that the cluster's radius R <C A <C r. Validity of that criterion permits one to use 
the local, dipolar approximation to describe the electrodynamic response of the cluster to 
the fluctuating electromagnetic fields that are responsible for vdw interactions. Based on 
this argument, the interaction between two clusters A and B may be written 

pO pO 

V AB (r) = -K^jpL (2) 

Here P A and P# are the static polarizability of cluster A and B, respectively, a subject of 
much attention in recent years0'0lllill. Now, one may ask a simple question: how different 
is this cluster-cluster interaction from a naive estimate, V es t, obtained from summing pairs 
of interactions between the constituent a, b atoms? The latter may be written as a product 
of the interatomic pair interaction and the number of such interactions between the two 
clusters: 

V est (r) = N A N B V ab {r) = -KN A N B ^- (3) 

We compute the ratio of the exact result Eq. (2) to this estimate and define the ratio as 
a "screening" function S AB : 

<T (r) - VAB{r) - P " a{Na) P " b{Nb) - f (N )f (N ) U) 
S AB (r) - ^ - ^ - f A (N A )f B (N B ) (4) 

This screening function is separable, a product of two independent screening functions 
f(N), each of which is the ratio of the cluster's polarizability to a nominal cluster polariz- 
ability, equal to the product of the atomic polarizability and the number of atoms. Because 
each of these / functions is less than one if N > 1, then S r AB < 1 if either of the clusters 
possesses more than one atom. Thus the true interaction is reduced relative to the naive 
estimate V est . 

Figure 2 presents the screening function for the case of a Na cluster interacting with 
a second Na cluster. The static polarizalility of the Na clusters P°(N) was taken from 



theoretical results, obtained using the jellium model and density-functional calculationsta. 
This simple model was adequate to explain the experimental polarizability data§. For the 
extreme case Na = 1 = N B , the screening function is unity, by definition. For increasing 
N, typically, the screening causes S r AB to decrease below one. Note that the reduction is by 
a factor near 0.4 for two 40 atom Na clusters and this is essentially the asympotic value@ 
for Na, Nb — > oo. The dependence on number is not monotonic, however, because of shell 
structure of the electronic states of the cluster. This nonmonotonic behavior makes the 
figure interesting; notice the shape of a "bat head" in the center of the Figure 2. 

We remark that the hypothetical case of clusters comprised of weakly polarizable, inert 
gases represents an instance for which the functions whose product determines S r AB are 
close to 1, i.e. P A = N A P®. Indeed, in the framework of local, continuum electrostatics, 
this relation is exact if the Clausius-Mossotti relation describes the static dielectric constant 
e of the material, since the dipole moment of a dielectric clusteril in an external field E 
is [(e — l)/(e + 2)]i? 3 E. Thus, the interaction between two such inert gas clusters is not 
significantly screened. 

B. Cluster-surface interaction 

Next, we address very briefly the case of a cluster (A) near the surface of a semi- infinite 
perfectly conducting metal, a simplified version of a more general adsorption problem. In 
that case, the cluster's interaction with the surface obeys a relation! 



Vcond = (5) 



z 



Here K' = 3hc/(16ir 2 ) and z is the separation between the cluster's center and the 
boundary of the surface . The ratio of this cluster-surface interaction to a naive estimate, 
based on Na individual constituent atoms of the cluster is analogous to the ratio S r AB for 
the cluster-cluster interaction: 

This behavior therefore involves the same function as the previous problem^!. 

Finally, we evaluate the retarded cluster-surface interaction in the case of a solid made 
of molecules of static polarizability and density n. We specialize the discussion to the 



situation where the product nP^ <C 1; this weak- screening limit is essentially the extreme 
opposite case from the perfect conductor considered above. In this new limit, the interaction 
V W eak can be evaluated by integrating the cluster-molecule interaction density —KnP^P^/r 7 
over the half-space z < in the presence of the cluster, centered at distance d above the 
surface: 

y^d) = -nKP^l i<oJ ^- f (7) 
The result of that integration is: 

v M = -JL Kn 3Ek = -™. h ™m (8 ) 

mt[ ' 10 d" 80jt d> { ' 

We may compute the ratio of this interaction to that of the perfect conductor with the 

same cluster (Eq. fl5|)): 

Vweak _ 23 ^.^pO 
Vcond 1 5 

Note that these interactions strengths are similar if nP^ « 0.2. Even for such a small 
polarizability, therefore, the interaction is comparable to that obtained with a perfect con- 
ductor. In the case of an Ar solid, in contrast, the ratio in Eq. @ is about 0.2. 

III. NONRETARDED INTERACTION 

We first recall that the nonretarded interaction between two atoms a and b separated by 
a distance r is given by0 

s~tab 

Mr) = (10) 

where the constant Cg b is 

'if, rco 

Cf = — duP a {iu) P b (iu) (11) 

7T JO 

Here, P a {oS) and Pb{oS) are the dynamic polarizabilities of the atoms, continued to imag- 
inary frequencies u = iu. The atomic polarizability has the form 

^M = -E^jS (12) 
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where fo n is the oscillator strength that measures the probability of the transition from the 
ground state to the excited state n, at frequency u 0n . 

A. Cluster-cluster interaction 

To compute the interaction between two clusters A and B, we use the same idea as in 
the last section of considering the clusters as large atoms (E < r < A). In that case, we 
may write the cluster-cluster interaction as: 



V, 



nr / \ 6 



C* B = — / duP A {iu) P B {iu) (13) 

7T JO 



We propose a simple form for the dynamical polarizability of the clusterii: 

P A (iu) = -Sjr (14) 

designed to give the static polarizability of the cluster at zero frequency, and the known 
asympotic behavior, Pa — > ZN "^e 2 '/ (mu 2 ) in the high frequency limit (free electrons) if we 
define uja with 

2 _ ZN A e 2 

Ua = ( } 

Here m is the electron mass and Z is the number of electrons in the atom. An analogous 
aproximation can be made for the atomic polarizability P a {iu) using a single characteristic 
frequency u a to determine its form. The characteristic frequency of the cluster, u>a, is larger 
than the corresponding result for a single atom, u a , by a factor of 1/^/Ja- We have tested the 
quality of our proposed function by comparing it with the dynamical polarizability of small 
metal clusters, calculated by Kresin within the random-phase approximation^. He finds two 
collective excitations (one surface mode and one volume mode) in terms of which he evaluates 
the dynamical polarizability. We have found that the overall frequency dependences of his 
expression and Eq. (|T4| ) agree very well, with a difference less than 0.5 %. 

Assuming the cluster polarizability to be that from Eq. ( |T^ ) (and using the correspond- 
ing approximation for the atomic polarizability), we compute the ratio S^ B between the 
nonretarded cluster-cluster interaction and the simple estimate that results from summing 
pair interactions, assuming Z — 1: 
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cw Cq B u a u b (u a + u b ) f— (P°)' 1/2 + (P b °)-V 2 

Figure 3 shows this function for the case of two interacting Na clusters. It shares the 
same qualitative behavior as the screening function in the retarded case (Eq. (|4])). In the 
case a = b, Na = Nb, S^b = f\ (Na)- This compares with the result f\ in the fully 
retarded limit (Eq. (f|)). 

B. Cluster-surface interaction 

The nonretarded vdw interaction between an atom and a surface at separation d is given 
by! 

^ya—surf 

V a . surf (r) = (17) 

with 

* f> r°° 

C a-surf = duP a {lu)g{lu) (18) 

47T JO 

Here, g(iu) is the dielectric response function of the substrate 



9{m) = 7uu)Ti (19) 



which can be approximated by 



9H = (20) 

if the response is dominated mainly by the resonance at u s (the surface plasmon resonance 
in the case of a metal). The parameter go is equal to 1 in the case of a free electron metal, 
in which case the ansatz Eq. fl2Dp is exact, with oj s = v2 times the bulk plasma frequency. 
Following the spirit of the previous sections, the van der Waals coefficient for the cluster- 
surface problem is 



Na 4tt Jo N a 8 (l + ^) 
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After carrying out the integrations for C 3 sur ^ and C3 sur ^ the screening function for 
the cluster-surface interaction results: 

^jA—surf _|_ Us 

SA-surf = ^a-surf = f A 1 7 u s A fl<= ( 22 ) 



We plot this function in Figure 4, for different ratios uj s /uja- As can be seen in the figure, 
SA-surf —> \^Ta, if Us 3> wa, and goes to when u; s <C uja- This relative interaction ratio 
is the smallest in the latter case (since J'a < !)■ 



IV. SUMMARY 



In this paper we have evaluated two kinds of interaction (cluster-cluster and cluster- 
surface) in two extreme opposite regimes of separation (nonretarded and fully retarded). 

We have described a particularly simple result for the fully retarded interactions involving 
clusters. The interaction between two clusters is reduced (relative to the naive estimates) 
by a screening function, derived from the product of individual cluster screening functions 
f{N), that depend on the electronic properties of the cluster. Because the latter shows 
interesting, nonmonotonic behavior as a function of atomic number, the dependence on the 
individual clusters' numbers is nontrivial, as exemplified in Figure 2. Qualitatively similar 
behavior occurs in the nonretarded case, seen in Figure 3. 

This work can be generalized in many ways. One potential extension involves explo- 
ration of the intermediate regime of separation. While straightforward in principle, this 
is complicated in practice, so we simply mention results obtained elsewhere for analogous 
problems in this intermediate separation regime (the dashed region in Figure 1). In the 
case of the interaction between two identical atoms, the effect of retardation is to reduce 
the nonretarded interaction by a factor of two when the atoms are separated by one-tenth 
of the characteristic wavelength in their excitation spectrum (i.e., r is about 40 nm if the 
relevant energy is 3 eV§). In the case of the interaction between an inert gas atom, or small 
molecule, and a surface, the corresponding distance is somewhat smaller, about 20 n m§. 
These distances are sufficiently small that one should not ignore the effects of retardation 
in many applications. If one were to accept the linear interpolation in Figure 1, the ratio of 
this interaction to the (extrapolated) nonretarded interaction at r = A would be 0.49. 
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The second problem of potential interest is the effect of a medium in which the clusters 
are dissolved. This is the principal focus of our future research. A third problem of potential 
interest is larger clusters (i.e. radii not much smaller than their separation). Both of these 
present no question of principle but neither problem is as simple computationally as the 
problems described here. 
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FIG. 1: Different regimes of the van der Waals interaction between two 20 Na atoms clusters, in 
units of a characteristic energy AE = Tiuq. The nonretarded regime (NR) when R/X <C r/X <C 1 
is analyzed in Section III, and the retarded regime for separations such that r/X S> 1 S> R/X 
is explored in Section II. The intermediate regime (dashed line) is not considered in this work, 
but note that it can be interpolated between the r -6 (nonretarded) and r~ 7 (retarded) functional 
forms. 
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FIG. 2: Retarded screening function S t (Na,Nb) for the case of two Na clusters containing 
and Nb atoms, respectively. 
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FIG. 3: Same as Fig. 1 for the nonretarded interaction. 
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FIG. 4: Nonretarded screening function SA-surf f° r the- case of a Na cluster of Na atoms interacting 
with a surface which has a characteristic energy Hoj s . 
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